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Abstract 



The strain energies in straight and bent single-walled carbon nanotubes 
(SWNTs) are calculated by taking account of the total energy of all the occu- 
pied band electrons. The obtained results are in good agreement with previous 
theoretical studies and experimental observations. The Young's modulus and 
the effective wall thickness of SWNT are obtained from the bending strain 
energies of SWNTs with various cross-sectional radii. The repulsion potential 
between ions contributes the main part of the Young's modulus of SWNT. 
The wall thickness of SWNT comes completely from the overlap of electronic 
orbits, and is approximately of the extension of tt orbit of carbon atom. Both 
the Young's modulus and the wall thickness are independent of the radius 
and the helicity of SWNT, and insensitive to the fitting parameters. The 
results show that continuum elasticity theory can serve well to describe the 
mechanical properties of SWNTs. 

PACS numbers: 61.48,+c, 63.20.Dj, 71.20.HK, 71.25.-s 
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Since their discovery in 19910, Carbon nanotubes (CNTs) have invoked considerable 
interest! in the last decade . There are many works on both the theoretical! - ! and the 
experimental! studies about the electronic structure of CNTs, and many exciting and novel 
properties have been discovered. For example, it was found that the insulating, semi- 
metallic, or metallic behavior depends upon the radius and the helicity of CNTsl On the 
thermal and the mechanical properties, the tubes are significantly stiffer than any material 
presently known@. To understand these many intriguing properties, many groups have cal- 
culated the strain energ yfi and the Young's modulusEJllB of single-wall carbon nanotubes 
(SWNTs). Among these calculations, many depend on the choice of an empirical potential 
between the carbon atoms, such as Tersoff- Brenner potential^. Lenosky et a/.0 employed 
an empirical model with three parameters reducible to a continuous model with two elas- 
tic moduliEl. They showed that the continuum elasticity model serves well to describe the 
deformation of multi-wall carbon nanotubes (MWNTs). Resent theoretical studies on the 
Young's moduli of SWNTs&ii show some discrepancies coming from the adoption of dif- 
ferent empirical potentials and different relations in the continuum elasticity theory (CET), 
especially, the different values of the effective wall thickness of SWNT. How to calculate the 
Young's modulus of SWNT is still an open question. 

Here we present a simple method for the computation of the strain energy of straight 
SWNTs directly from the electronic band structure without introducing any empirical po- 
tential. This method had also been extended to calculate the strain energy of bent tubes. It 
is found that the wall thickness of SWNTs can be calculated simply from the band electrons, 
and the Young's modulus by consideration of both the repulsion energy between ions and 
the bond-length dependencies of the electronic energy. Our results show that CET can well 
describe the bending of SWNTs and that both the Young's modulus and the effective wall 
thickness are independent of the radius and the helicity of the tubes, and insensitive to the 
fitting parameters. We obtained the Young's modulus of SWNT about 5 TPa, 5 times larger 
than the value of MWNT or graphite bulk samples, and the effective wall thickness about 
0.7 A, the size of carbon atom. 

Generally, the total energy of the carbon system is given by the sumElil!: 

E total = E e i + E rep , (1) 

where E e i is the sum of the energy of band electrons of the occupied states and E rep is given 
by a repulsive pair potential depending only on the distance between two carbon atoms. 
They are given by 

E e i = E k , (2) 

occ 

and 

E rep = ^2^24>(r ij ), (3) 

i j>i 

, respectively. Since 0(r) is a short-range potential^, only interaction between neighbor 
atoms needs to be considered. On account of the relaxation effectHi, the bond length of 
SWNT is slightly larger than that of graphite (r = 1.42 A). However, even in Cqq, for 
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which the relaxation effect is significant on account of its small radius, calculations show 
that the energy contribution of the bond-relaxation can still be safely ignored^. The total 
energy can now be rewritten as: 

where the first term on the right hand side of E to tai is the sum of the repulsion energy 
between ions and the electronic energy contribution of the bond length change with br^ 
as the change of the distance between the zth and the j'th atoms in SWNT from that in 
graphene. The second term is the electronic energy contribution of the angular change of the 
bond, when rolling from graphene to SWNT. The positions of the atoms of straight SWNTs 
are located on the cylindrical surface of the tube when the relaxation effect of the bonds is 
neglected, br^ is proportional to p~ 2 , where p is the cross-sectional radius of SWNT, and 
the first term of Eq. (f|) can be ignored, since it is of p~ 4 order. Therefore, the strain energy 
of straight SWNTs comes from the curvature-induced electronic energy change, and can be 
obtained by taking account of the electronic energy of all the occupied bands. 

In order to calculate the electronic energy bands of SWNTs, we use a simple nearest- 
neighbor tight-binding (TB) model. This model contains nine TB parameters of graphite: 
Four hopping including, V ssa = —6.679, V spa = —5.580, V ppa = 5.037, V pp7T = —3.033 in unit 
of eV; four overlapping integration including, S SS(T = 0.212, S spcT = 0.102, S ppa = —0.146, 
Sppn = 0.129 and an energy difference between the 2s orbit and the 2p orbit of the carbon 
atoms A£=(£ 2s ~^2 P ) = —8.868 eV.i The model has been used widely for the calculation of 
the electronic properties of both graphenes and SWNTs. In general, these TB parameters 
depend upon the bond-length in the way@, 

v xx'^( r ) = v xx'^( r o) ■ exp(- 7 (r - r )). (5) 

However, in case of straight SWNTs, the p~ 2 order dependence of the strain energy will not 
be affected even if we ignore simultaneously these dependencies and the repulsion energy. 

With the notation used by White et a/.i, each SWNT is indexed by a pair of integers 
(ni, n 2 ) corresponding to the lattice vector R=n\di + n 2 a 2 on the graphene, where a 1; a 2 are 
the unit cell vectors of the graphene. The tube structure is obtained by a rotation operation 
Cn and a screw operation S(h,a). The operation Cn is a rotation of t£ about the axis, 
where N is the largest common factor of n\ and n 2 . The S(h, a) operation is a rotation of 
an angle a about the axis of SWNT in conjunction with a translation of h units along the 
axis, which both h and a depending on the tube parameters!. Let [m, I] denote a primitive 
unit cell in the tube generated by mapping the [0, 0] cell to the surface of the cylinder first 
and then translating and rotating this cell by I applications of the rotational operator Cn 
followed by m applications of S(h,a). Because S(h,a) and Cn commute with each other, 
we can generalize the Bloch sums, and obtain the Hamiltonian matrix: 

H$ A (k,n)=H% B (k,n)=£ i 5 ij , 
Hf 1 A (k,n) = (Hf(k,n))\ 

Hf(k,n) = ^^[yAlW+ifcMr)]^), (6) 
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where [Am(r), Al(r)] (r = 1,2,3) are the cell indices of the primitive unit cells located by 
the three nearest-neighbor atoms B of atom A in the tube, A and B are two independent 
carbon atoms in a primitive unit cell of SWNT. Let n = 0, 1, . . . , N — 1 represent the N 
sub-Brillouin Zones, A; be a one-dimensional wave vector. E\ = £ 2s , £% = £ 2p for i = 2,3,4. 
Taking the 2p wave function as a vector and 2s wave function a scalar, one can easily obtain: 

^ P] {r)=(e Ai -e B]i:r) )V PTm 

-(eAi ■ U r )(e B .(r) ■ u r )(V ppn - V ppa ), 
V**(r)=(e Bi (r)-u r )V spa , 

■ Ur)V spa , 

V s A f(r)= V ssa . (7) 

where ii r (r = 1, 2, 3) be the unit vector from the atom A to its three neighboring atoms B. 

and e B .(r) are the unit vector of the 2pi wave function of atom A and the unit vector of 
the 2pj wave function of atom B, respectively. The overlapping integration matrix S(k,n) 
has the same form as the Hamiltonian matrix, with four overlapping integration parameters 
to replace the four hopping integration parameters, and with unit to replace the energy 
of 2s and 2p wave function £ 2s and £ 2p - Thus we obtain an 8 x 8 Hamiltonian matrix 
7i(k,n) and an overlapping integration matrix S(k,n). By solving the secular equation 
Ti(k, n)Ci(k, n)=EiS(k, n)Ci(k, n), we can calculate the electronic energy band Ei of SWNTs. 

Taking account of the total energy of all the occupied band electrons in SWNTs relative 
to that in graphene, we have calculated the strain energy E s of the straight SWNTs. With 
the possible bond length dependence of the TB parameters being neglected, and with the 
real bond length of SWNT \u r \, where u r represents vectors between the nearest neighbor 
atoms in the tube@, we have calculated the direction cosines • u r of Eq. (|7|). Fig. 1 shows 
that E a depends only on the radius p of the tubes. The characteristic behavior E s =C/p 2 is 
found with C ~ 1.44 eVA 2 /atom, in good agreement with previous calculated value 1.340 
or 1.53 eVA 2 /atorrJll, and excellently close to the value of 1.57 extracted from the measured 
phonon spectrum of graphitelll. 

Recently, "curved SWNTs" and "torus-like SWNTs" have been foundS They still have 
the sp 2 bond structure, but it is predicted to have pentagon-heptagon defects^. In "curved 
SWNTs", the bond-length is nearly the same as that in graphite sheet, since the distortion 
that is created by the bending nature of the curved tube is topologically relaxed by the 
inclusion of fivefold and sevenfold rings. However, the application of an external force 
moment at the two ends of the tube gives a different deformation. The hexagonal structure of 
the tube will not change until it reaches a critical bending curvature!!. The tube undergoes 
only a simple compression on the inner side, and a stretching on the outer side. In the 
following discussion, the word "curved SWNT" refers to the growthed SWNT with pentagon- 
heptagon defects, and the word "bent SWNT" refers to SWNT that bends with outer- 
stretching and inner-compressing deformations under external force moments applied to the 
two ends of the tube. Using an empirical model employed by Lenosky et a/JHl, Ou-Yang 
et alM have developed a macroscopic continuous elastic model to calculate for "curved 
SWNT" . In their work, the strain energy of "curved SWNT" come from the angular change 
of the bonds or the curvature of the tubes. However, in the case of "bent SWNT" , the bond 
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length effect will contribute the main part of the strain energy. In what follows, we will 
treat only the latter case. 

The "bent SWNT" surface can be described byEl 

Y(s, 0) = f(s) + p[N(s) cos + b(s) sin 0] (8) 

where f(s) is the position vector of the axis, and < s < I, is the arc-length parameter along 
the bent SWNT axis. < < 2tt. N(s) and b(s) are unit normal and unit binormal vector 
of r(s), respectively. The position of each carbon atom is described by two parameters, s 
and 0. The two operations Cn and S(h,a) can still be used to determine the positions of 
atoms in the SWNT. Therefore, a translation of h units along the axis of SWNT means 
an addition of h to s, and a rotation of a about the axis means an addition of a to 0. 
However, because the rotational symmetry about the axis of the bent SWNT is broken, Cn 
and S(h,a) are not symmetry operations. It is necessary to generalize the Bloch sums in 
the crystal unit cell containing M. x iV primitive unit cells of SWNT. Here M. is the length 
of the cell along the axis direction( unit is h ). M. = 1{n 2 + n^ 2 + 11^2) /N 2 , for ni — no 
not a multiple of 3iV, and M = 2{n 2 + n 2 2 + n!n 2 )/(3iV 2 ), for n x — n 2 a multiple of ZNti. 
We calculate only SWNT with constant radius of curvature R. It is not principal difficult 
to extend the present treatment to general bent SWNTs. In bent SWNT, V^ B (r) depends 
on the position of the atom A, and will be written as V^ B (l,m;r), where [m, I] are indices 

of the primitive unit cell of the atom A. When SWNT is bent to a different direction N(s), 
V{j B {r) will be different, but the Hamiltonian matrix elements are almost independent of 
the bending direction. We have found that the anisotropy effect is very small. Similar to 
the case of the straight SWNT, it is easy to obtain the (8 • M. ■ N) x (8 • M. ■ N) matrices 
of the Hamiltonian and the overlapping integration. 

Since the change of bond length 5r in bent SWNT is proportional to the energy 
contribution of the bond stretching and bond compressing will be of the order of It is 
necessary to calculate both the electronic energy E e i and the repulsion energy E rep between 
the ions. In order to fit the force constant of graphite@, we take 7 = 1.024 A" 1 , 0' = ^ = 

—13.637 eV/A and " = 60.4 eV/A 2 . With these parameters, we arrive correctly the second 
derivative of the stretching energy E c of SWNTs be D = ^ = 58.5 eV and the Poisson 
ratio be o = 0.24, where e is the relative compression along the axis of SWNT. 

Fig. 2(a) shows the strain energy Eb per atom of the (5,5) SWNT as a function of the 
bending radius R. The data follow quite well with the expected behavior Ej, = E s + X/R 2 . 
A least square fit to the data yields a value of A ~ 173 eVA /atom. Previous studies on 
"curved SWNTs" give a simple formula can be givenEl, 

CR C C . . 

b = J7W = T 2 ^7 + 2-R 2 - (J) 

In comparing with our results, we find that the value of A for "curved SWNT" is equal to 
C/2, only 0.7 eVA /atom. It implies that the strain energy of pentagon-heptagon defects is 
far less than the strain energy of the stretching and compressing of the bond length. The 
experimental fact that the deformations of "bent SWNT" are change of bond length rather 
than the pentagon-heptagon defect reveals that there is a high potential barrier between the 
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two deformations to prevent the change of the hexagonal structure under the addition of a 
moment of external force at the two end of SWNT. 

By CET, we can calculate the Young's modulus Y of SWNT from three different strain 
energies. They are the rolling energy E s , the compressing or stretching energy E c and the 
bending strain energy AEf,. The three energies are given by 

E s = (10) 
P 

E c = \De\ (11) 

A£ 6 = A (12) 

where e is the relative stretch or compression along the axis of SWNTs, and D is the second 
derivative of E c . The three quantities C, D and A are given by 

c= w L ^ Y - b3 ' (13) 

D = n-Y-b, (14) 
A = fr-6(p 2 + 6 2 /4), (15) 

respectively. Where Q = 2.62 A /atom is the occupied area per carbon atom in SWNTs, b 
is the effective wall thickness of SWNT. Previous calculation indicate that the value 

of D of SWNTs is about 58 eV/atom, same as that of graphites. However, since the wall 
thickness is not well defined in single-layered structure, various values of b are used in the 
studies, thus the obtained Y are quite different. Lu@ and Hernandez et a/.0 took the 
interwall distance of graphite (3.4 A) as the thickness, and obtained the average Young's 
modulus of SWNT about 1 TPa, in consistency with the corresponding measurement in 
multiwall nanotubesll and bulk graphite samples. But the average value of Y can not 
describe all kinds of deformations of SWNTs, such as the rolling of graphene and the bending 
of SWNT, though it can describe the stretching and compressing deformation along the axis 
direction of SWNT. Yakobson et a/.0 have given Y = 5.5 TPa and b = 0.66 A by using 
the rolling energy formula of graphite sheet [ Eq. (0)] and the stretching energy formula 
of graphene or SWNT [Eq. (|13"D 1, simultaneously. The obtained value of b is about the n 
orbital extension of carbon atom, which corresponds to the general fact that elasticity results 
from the overlapping of electron cloud between atoms. However, since Eq. ( |T3"D describe 
the rolling of single- layered graphene, the results given by Yakobson seems to correspond to 
graphite sheet rather than SWNT. The Young's modulus and the effective wall thickness of 
SWNT and their dependence of the tube radius and helicity still remain unknown. 

With Eq. (|15D , from the calculation of the bending strain energy of SWNTs with various 
radius, one may simultaneously find Y and b of SWNTs. Fig. 2(b) shows the relationship 
between A and the radius of the tube in the form A = b*p 2 + a*. In comparison with the Eq. 



(15), it implies that both Y and b are independent of the radius and helicity of SWNTs. The 
value of b* = 15.3 eV/atom is consistent with the value D/4 = 14.7 eV/atom. However, it 
is difficult to obtain the exact value of a*, because the first term of A is much greater than 
the second term a*, to introduce high errors in the our fitting of a*. Carefully analysis of 
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the strain energy of "bent SWNT", gives only that the electronic energy from the angular 
change of the bond can contribute to the p° order of the bending strain energy of SWNT. 
The other terms, including the repulsion energy between ions and the electronic energy from 
the nonzero 7 effects (Eq. depend on Sr, hence on the radius of the tube p. For wall 

thickness, is is unnecessary to consider the bond length dependence of TB parameters and 
the repulsion energy. It require only to calculate 7 = bond angular contribution E e io of 
the electronic energy. When 7 = 0, the -^th order perturbation of Hamiltonian is zero, and 
E e io = The p° order term of A come completely from the A e /o and the residual part of 
A affects only the value of p 2 order. The exact value of a* can be obtained by calculating 
E e i . Fig. 2 (c) shows the expected relationship E e i = X e i /R 2 . Fig. 2(d) gives values of 

A e ; for some SWNTs. It lead to a* = 1.05r 2 eVA 2 /atom by fitting it to a* + a%p + a 2 p 2 .@ 
The wall thickness of the tube is supposed to be identical to that of graphite, then from Eq. 
(0) and Eq. QTJ), it gives a* = |(1 - o 2 )C » 1.0r 2 eVA 2 /atom. Therefore, b is about 0.74 
A, and Y is about 5.1 TPa. This shows that both Y and b are independent of the radius 
and helicity of the tube, and the Young's modulus of SWNT is five times greater than the 
average value of MWNT. The obtained value of b is independent of the fitting parameters 
7 and , and Y is also insensitive to these parameters. 

In summary, our calculation shows the following results: the strain energy of the straight 
SWNT come mainly from the occupied bands electrons, The obtained Young's modulus of 
SWNT is independent of the radius and the helicity and is much larger than the modulus 
of the bulk sample. The effective thickness of SWNT is about the size of the carbon atom, 
far less than the distance between the layers of the graphite. These results show that CET 
can well describe the deformation of the bent tubes. 

The authors acknowledge the useful discussions in our group. We would like to thank 
Dr. X.-J. Bi, Mr. Y.-H. Su and G.-R. Jin for correcting an earlier version of the manuscript 
and thank Prof. Y.-Z. Xie and Dr. H.-J. Zhou for correcting the manuscript in English. 
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FIGURES 
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FIG. 1. The strain energy per atom versus the radius of (n, m) tubes, where n = 6 ~ 13 
and m = 0,1,2 or n. The solid line corresponds to a least square fit to the ^ behaviors. 

C = 0.71 • ro 2 ~ 1.44 eVA 2 /atom. C is independent on the helicity of tubes. The p~@ behav- 
iors is clearly shown in the inset ((3 pa 2.03 ). Here ro is 1.42 A, the cohesive energy of graphite is 
-7.39 eV. 
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FIG. 2. (a) Strain energy per atom versus the bending radius R in (5,5) tube. The solid line 
is a fit to the E s + X/R 2 , where E s is the cohesive energy of straight (5,5) tube . A = 86.1 x 1.42 2 
eVA /atom, (b) the value of A of some (n,0) and (n,n) tubes. The solid line is a fit to a* + b*p 2 , 
b* ~ 15.27 eV. (c) E e iQ versus the bending radius R in (5,5) tube. The solid line is a fit to 

+ X e io/R 2 . The zero point of E e i has moved, (d) The value of A e ;o of these tubes. The solid 
line is a fit to ao + a±p + a2p 2 , and the ao ~ 1.05. 
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